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ABSTRACT 


The viewpoint is taken that in large portions of the atmosphere, wind speeds (V,,), divergence (D,,), and vor- 
ticity (¢,,) obtained under the gradient wind assumption are considerably more accurate representations of true 


conditions than those one obtains from the geostrophic assumption. 
The equation for the vorticity of gradient winds has a considerable resemblance to “balance’’ equa- 
Gradient winds and their space derivatives may have advantages over other 
The quantities V,,, D,,, and ¢,, may be computed quite easily from 


204 f er. 
tions of Charney [3] and Phillips [17]. 
formulations due to their relative simplicity. 


geopotential data at grid points by use of high-speed computers. 


and theoretical meteorology are suggested, 


1. INTRODUCTION 


The geostrophic wind equation, an empirical relation 
which equates horizontal pressure and Coriolis forces, is a 
well-known cornerstone of both theoretical and applied 
meteorology in extratropical latitudes. One might expect 
that a generalization of the geostrophic equation to 
include a third term of major importance, namely, the 
centrifugal force, would lead to more accurate modeling 
of winds than the simple geostrophic relation provides. 
However, the reader will immediately realize that the 
centrifugal force V?A, where V is wind speed and K is the 
trajectory curvature, is small in regions of weak winds, 
and that in regions of straight flow the gradient wind 
automatically degenerates into the geostrophic. For 
example, Neiburger et al. [13] compared geostrophic and 
gradient winds with winds on two 700-mb. 
charts having a mean observed speed of 18 kt. and found 


observed 


that gradient winds were not significantly better than 
geostrophic. In a study of strong winds (having an 
average speed of 94 kt.), Endlich and McLean [7! found 


Equations are derived for computing V,,, 


Possible applications of gradient winds in practical 


that in cyclonically-curved jet streams, geostrophic winds 
overestimated wind speeds measured by aircraft by 28 kt. 
on the average (18 percent of the geostrophic speeds of the 
sample) while gradient winds differed only minutely 
In other words, the centrifugal 

The standard 


from measured values. 
force was 18 percent of the pressure force. 
deviation of geostrophic departures (after accounting for 
observational errors) was estimated as 22 kt. while the 
standard deviation of gradient departures was estimated 
as 5.5 kt. From these data, one would conclude that in 
and around jet streams, and probably also in relatively 
strong currents at other altitudes, consideration of the 
centrifugal force accounts for the major part of the 
ageostrophic components. Remaining ageostrophic com- 
ponents are due to cross-contour flow, to vertical motions, 
and to friction. However, studies of constant-level 
balloon trajectories (Neiburger and Angell [12], Angell 
[1]) indicate that cross-contour components are relatively 
large, being on the average about two-thirds of the magni- 
tude of the along-contour (centrifugal) components. 
Thus, the question of the relative magnitude of the two 
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ageostropliic components is presently unresolved. Never- 
theless, Angell estimates that the percentage error in 
estimating winds geostrophically exceeds 29 percent half 
f the time, while the error of the gradient assumption ts 
‘considerably less, exceeding 11 percent half of the time. 
In practice, one can include the centrifugal components 
by use of gradient winds; however, the cross-contour com- 
ponents, though important in initiating disturbances and 
in subsynoptic phenomena, are very difficult to estimate 
from conventional data. 

The hypothesis is therefore proposed that gradient 
winds provide an estimate of actual winds which is 
sufficiently accurate to make them useful in describing 
the dynamies of synoptic-scale processes. In support 
of this hypothesis, one may recall that the propagation 
ol evelone Waves Wis explained (see, for example, the 
review by Bjerknes [2]) on the basis of divergence of 
gradient winds in advance of upper-level troughs and 
Moreover, Phillips [17] has shown 
“meteor- 


convergence behind. 
that strictly 
ological noise” and are not suitable for computing initial 


non-divergent winds introduce 
data for numerical prediction based on integration of the 
primitive equations. Gradient winds, which are divergent 
vet approximately geostrophic, might provide a simple, 
relatively noiseless balance condition between winds and 
the geopotential field. It will be shown later that the 
equation for gradient vorticity (€,,) Is quite similar to 
the well-known balance equation (Petterssen {14}, Charney 
[3)) and to a vorticity equation based on a divergent 
wind (Edelmann [6], Phillips [17]). 

One of the main difficulties to be expected in computing 
winds (and especially divergence fields) from present 
upper-air data is that random errors caused by observa- 
contained the results. 


tional inaccuracies will be 


Therefore, smoothing of computed values may be a 
prerequisite to their use. 

In order to determine whether advantages accrue, in 
practice, from the explicit calculation of gradient) winds 
and their divergence and vorticity fields, in the future it 
will be 
parison with observed winds and with divergence and 
In the present 


necessary to carry out computations for com- 


vorticity calculated by other techniques. 
article, a simple method of successive approximations ts 
developed for computing V,,, Dy, and ¢¢,. The nature 
of the computational scheme is such that it can be carried 


out quite simply and rapidly by use of electronic computers, 


2. COMPUTATION OF VELOCITY, DIVERGENCE, 
AND VORTICITY OF GRADIENT FLOW 


The basic assumptions of gradient flow (e.g., Holmboe 
et al. [11]) are that air parcels maintain constant speed 


(i.e., do not cross the contours on constant pressure 


charts toward higher or lower heights) and move hori- 


zontally and without frietion. The equation which 


results on introduction of these assumptions into the 
horizontal equation of motion may be written as follows: 
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(d9/ds) n=—f(V,—V,) n 1) 


where V, is the gradient speed; @ is the wind direction 
taken as positive in the counter-clockwise direction; s 
is distance along the path; d@ds is AK the curvature of 
the path: n Vo Vo 
left of the path, V operates at constant pressure, and 


is the geo- 


is a unit vector directed to the 


o=g2; f is the Coriolis parameter; and V, 
strophie wind speed. The term on the left side represents 
the centrifugal force which is proportional to the differ- 
and V’,. 


If we post-multiply equation (1) veetorially by the 


ence between V,, 


vertical unit vector k we obtain 


IN. IV.- (19 ds)(n x k). (2) 
Also 3) 
nxk=(1/ Vo 


Therefore V,,=|1—(V,,7/ Ve )(da/ds) (4) 
After a further simplification, this form of the gradient 
wind equation will be used to compute V,,.. The simplifi- 
cation involves the use of Blaton’s formula 


dé dt da (= 


os 


a aa 


Oop) V os 


i.e., We neglect the term containing the vertical advection 
of 
sional trajectories tend to be less evclonic than trajectories 
The differ- 
ence is greatest in fronts or other phenomena where 
It is of lesser, 


If vertical motion is present, actual three-dimen- 
estimated on isobaric charts (Danielsen [5]). 


vertical motions reach maximum values. 
and perhaps negligible, importance in larger-scale motions. 
As there is no simple method for estimating the magnitude 
of the difference, it will be neglected in the remainder of 
this paper. 

We will now define a parameter k= (1) )(déds) or 
from equation (5) 


k=(1/ Vo )(09/ds) + (1/ Ve (6) 
where & represents the trajectory curvature, the first term 
on the right measures the contour curvature, and_ the 
second term measures the local change in wind direction. 
For cyclonic trajectories k>O and for antieyclonic, 
k<0. Since we have assumed flow parallel to contours, 
04 Of may be evaluated as the local change in contour 
direction, This direction may be measured counter- 
clockwise from east so that 


(r,/u,)=tan~! (—¢,/¢,) 7) 
where the subseripts # and y denote partial derivatives 
eastward and northward. 

To evaluate equation (6) we use 0@/0s=cos 6(08 dr) 
+sin 6(00/0y) where @ is given by equation (7) and, as 


a matter of convenience, substitute VW, for VV. Then 
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(o,” | 26,bybr,4 odyy) — Ort ) |. 
(S) 


Given the distribution of @ at grid points on a constant 
pressure chart, one can easily compute a finite difference 
approximation to the first term in the bracket. A short- 
ened version of this term has been used by Phillips [16| 
in computing contour curvature. The second term in the 
bracket can be evaluated as the difference between geo- 
potential gradients on charts separated in time, or from 
the space variations of geopotential tendency. 

For progressive, wave-shaped contour patterns, the 
curvatures of trajectories and contours are of the same 
sign (unless the waves are traveling faster than the winds) ; 
however, the trajectories are less sharply curved than the 
contours.' In other words, the term (1/17)(06/0f) is 
generally of opposite sign but smaller magnitude than 
0# 0s. The smoothing effect of a finite-difference evalua- 
tion of 080s tends to reduce its magnitude toward the 
magnitude of the trajectory curvature. Therefore, on 
upper-level charts, especially where V is large, it may 
prove feasible to neglect the term containing (1/V)(0@ df). 

Having obtained a value of & at each grid point by use 
ol equation (S), We may write equation (4) in vector or 
scalar form as 


V.,=(1—kV,,2)V, or V,,=(1—kV,,2) Vy. (9) 


This quadratic equation mas be solved for V os since k 
and Vy are known. In the past, the solution has ordi- 
narily been obtained graphically; e.g., Gustafson [9]. If 
the equation is solved by electronic computer, the square 
root in the quadratic formula is calculated by use of an 
initial guess and successive approximations. It appears 
that a solution to equation (9) may be obtained more 
directly and simply by the following process of successive 
approximations. As a first guess, we may substitute 
V- for V,2 on the right side of (9). Then if we replace 
the subscript on V7), by the number of the approximation 
we obtain 


Vi=(1-kV2)V, 


(10) 


The convergence of V, to Vy, is quite rapid as shown 
in fizure 1. Here we have assumed a typical jet stream 
condition where V,=125 kt. and k=2x107° kt.-* (ey- 
clonic curvature) so that V,,=100 kt. The approxima- 
tions fall alternately on opposite sides of V,, and in this 
particular case, VW; approaches V,, within 1 percent. For 


relationships between streamlines and trajectories are discussed by Holmboe 
1} and Haltiner and Martin [10], among others 
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Figure |.—Graph showing the convergence of successive approxi- 


mations to true gradient wind speed (V,,) using the geostrophie 
speed (V,) as the initial guess. Solid lines join approximation 
for cyclonic curvature; dashed curve joins approximations for 
antieyelonic curvature, Approximations computed using equa- 
tion (10) in text 


V,=83.3 and & 107° (antievelonic curvature), the 
approximations approach V,, asymptotically. this 
(100 kt.) within | percent. 
It can be shown that for these two examples, convergence 


example, Vy approached V, 


will occur if the absolute value of the initial guess is less 
than 500 kt. In general, the region of convergence is 
quite broad. The geostrophic speed is a convenient 
initial guess which is in the neighborhood of the gradient 
speed and therefore results in an acceptable estimate of 
the gradient speed after a small number of successive 
approximations. 

For anticyclonic flow and for a given radius of curvature, 
up to a limit of 2V,. Beyond this 
point, the quantity 44V<—1 and equation (9) has 


increases with V, 
complex solutions which have no physical reality. Gra- 
dient wind balance does not exist in these “unstable 
ridges” and the approximations of equation (10) will not 
converge. Grid points at which this equation diverges 
may be marked with a special symbol to indicate the 
absence of gradient flow. Unstable ridges, though rather 
uncommon, are synoptically important since cut-off 
Lows sometimes form downstream (Bjerknes [2}). 

We will assume that by using equation (10) we have 
obtained a satisfactory estimate of V,, at each grid point 
(with the exceptions mentioned above). Now it. is 
convenient to define a quantity 


B=—-kV,,/ 
so that equation (9) becomes 
V,.=(1+B)V, or (11) 


B is a non-dimensional term which expresses the departure 
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of the gradient speed from the geostrophic as a fraction 
of the geostrophic speed. It may also be thought of as 
the correction to be made to the geostrophic wind in 
order to obtain the gradient wind. For convenience, 
it will be referred to as the “gradient correction.”” In a 
roughly sinusoidal pattern of contours in which a jet stream 
is imbedded, one would expect centers of negative B 
with central values of about 20 percent where the jet 
crosses troughs, and similar centers of positive B where 
the jet crosses ridges. 


The “oradient divergence” is defined as 


BV -NV,. (12) 
The divergence therefore depends on the geostrophie wind 
and the field of the gradient correction B. To obtain a 
quantitative estimate of the first term on the right side 
of equation (12), we take a typical value of geostrophic 
wind speed as 30 m. sece.-' To obtain a typical value 
of VB, we assume that B changes from —0.20 in a trough 
to zero at an inflection point over a distance of 600 km. 
Then V,-VB equals 107° sec“! The second term on the 
right side of (12) is generally an order of magnitude 
smaller but may be of importance in strong meridional 
currents in middle or low latitudes. Therefore D),,~ 107° 
sec.~' which is of the order expected for progressive waves. 

The relative vorticity of the gradient wind at any 
point is 


V,=1+ TBXV, (13) 
where is the geostrophic vorticity. Since Bis 
negative in cyclonic flow and positive in anticyclonic, the 
term (1+ 8) on the right side of (13) tends to make ¢,, 
either less or more than ¢,.  Petterssen [14] has shown 
that the geostrophic vorticity overestimates the vorticity 
There- 
fore, B provides a correction in the proper direction. 

of oy and equation (3), 
On expansion, the last 


of eyclones and underestimates it in anticyclones. 


term in (13) is 
Using these relations, equation (13) becomes 


Letting £6 using 
(V*o+ Bu,). 
v,B,—u,By. 

(1+ 1+ B)pu,+fr,B,—fu,B,=0. (14) 

In a machine computation of gradient winds and their 
divergence and vorticity, one might proceed as follows: 

1. Tabulate the geopotential at grid points on a 
constant pressure chart. 

2. Using geopotential data at neighboring grid points 
as finite-difference approximations to partial derivatives, 
compute w,, v, and V, (equation (3)), and @ (equation (7)). 

3. The curvature parameter & (equation (8)) may be 
Compute the first term in the 


computed in two steps. 
bracket from the grid point data of step 2. The second 
term in the bracket may be computed as the difference 
between past and present values of @, and @, or as the 
eastward and northward changes in geopotential tendeney. 
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4. At each grid point, compute V,, from V, and & 
(equation (10)). 

5. Compute B at each point from V,, and V, (equation 

6. Approximate VB by finite differences and compute 
D,, at each point (equation (12)). If desired, vertical 
motions may be obtained from J),, by use of the continuity 
equation with appropriate boundary conditions. 

7. Compute ¢,, from equation (13) or (14). 

It should be noted that these calculations do not require 
solution of first or second order operators over the grid. 
A rough estimate indicates that all of the above quantities 
(except vertical motion) can be calculated by a high-speed 
computer such as the IBM 704, for a 2000-point grid in 
approximately 30 seconds. 


3. COMPARISON OF THE VORTICITY OF GRADIENT 
WINDS WITH OTHER BALANCE RELATIONS 


Equation (14) bears a strong resemblance in its first 
three terms to corresponding terms in the “balance” 
equation (Petterssen [14], Charney [3]) written as 


put uZ+2r,u,+r7=0 (15a) 
or in terms of a streamfunction as 
— BY, — 2 Ver") 0. (15b) 


A slightly different equation has been obtained by 
Edelmann [6] and Phillips [17]: 


Bu, —2f (). (16) 


It should be noted that the vorticity specified by this 
equation generally differs somewhat from that given 
by equation (15). 

The similarities of equations (14), (15), and (16) are 
evidently due to the fact that they are all obtained from 
the horizontal equation of motion by approximating the 
acceleration in one way or another. The gradient wind 
approximates dV dt by the centrifugal force and is equiva- 
lent to expressing the geostrophic departure V' as 

V;,,? (da/ds). (17) 
As mentioned earlier, equation (14) is simply the vorticity 
of the gradient wind. Equation (15) is obtained by taking 
the divergence of the horizontal equation of motion and 
then neglecting the terms d(V - VW)/dt and Ve - (OW/dp). 
The wind field specified by this equation is non-geostrophic 
and non-divergent. Equation (16) gives the vorticity of 
a divergent, quasi-geostrophic wind which is obtained by 
expressing the geostrophic departure VW’ as k (dW/d/) 
and by approximating dV/dt geostrophically (Philipps 
[15]) so that 


V’=f-'k& (OV,/dt+V, + TV,). (1s) 
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A- a consequence of his demonstration that non-divergent 
intial wind fields introduce “meteorological noise” into 
iniegrations of the primitive equations, Phillips [17] has 
proposed obtaining initial wind data by solution of three 
siiultaneous equations (not repeated here) which give D, 
the horizontal divergence. ¢ is determined from equation 
and and with boundary conditions determine a 
wind field which is relatively “noiseless.” 

In the writer's opinion, the relative merits of expressing 
V’ by equation (17) or by equation (18) are difficult to 
judge on theoretical grounds. It follows that equations 
(14) and (16) are equally plausible; however, gradient 
winds and their divergence and vorticity appear much 
simpler to compute than the comparable quasi-geostrophic, 
divergent quantities. It may be logical, therefore, to 
compute fy, Vy, and ¢ by using equations (14), (15), 
and (16), respectively, from actual data for a number of 
cases. The three vorticity (or corresponding wind) fields 
might then be compared and tested in actual forecasts. 


4. APPLICATIONS OF GRADIENT WINDS 


The following applications of gradient winds appear to 
warrant further investigation: 

|. Specifving winds on current or forecast charts. 
Gradient winds might be used in flight forecasting for 
aircraft, in determining trajectories, and in locating jet 
streams (Endlich and McLean [8]) and fronts. 

2. Use as a balance condition between winds and the 
geopotential field; for example, computing initial 
vorticity fields for use in forecasting with the vorticity 
equation or in computing initial wind data for use in 
integrations of the primitive equations. 

3. Determination of a first approximation to horizontal 
divergence fields and corresponding vertical motions. 

$. As a means of introducing wind data into the objec- 
tive analysis of fields of geopotential or streamfunetion in 
place of the current “geostrophic” use of observed winds 
Cressman [4]). A) possible method of introducing ob- 
served winds according to the gradient assumption would 
be to compute the gradient correction (2) that corresponds 
to the geopotential field used as the initial guess. Ob- 
served winds divided by (1+ 8) would give geostrophic 
witids for use in the objective analysis procedure. 

As mentioned in the Introduction, the application of 
gradient winds to the problems mentioned above is based 
on the empirically supported hypothesis that they provide 
a considerably more accurate approximation to true 
conditions than the more restrictive geostrophic assump- 
tion provides; however, this hypothesis must be tested 
further with observational data before a final judgment 


of its validity can be made. 
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ABSTRACT 


\ practicable solution of the Penman formula for potential evapotranspiration by use of graphs is presented. 


For ease in computation, three graphs are used, 


These graphs, designed for Columbia, 8.C., can be adapted for 


any location by relabeling or correcting two of the families of curves. 


1. INTRODUCTION 


Potential evapotranspiration is a measure of the maxi- 
mum possible water loss from an area under a specified 
set of weather conditions. This phenomenon may be 
defined as the water loss from a vegetated surface which 
Penman {1} 


has shown that by the simultaneous solution of a Dalton- 


is supplied with adequate water at all times. 


tvpe equation and the energy-balance equation one can 
obtain a fairly simple, vet reasonably sound, equation for 
The general 


estimating potential evapotranspiration. 


form of the Penman equation is 

: 

D+F 

where /, is potential evapotranspiration; J) is slope of the 
saturation vapor pressure vs. temperature curve at the air 
temperature; Fis a constant (we use 0.27) in the equation 
R=F (T,—T,) where the Bowen ratio; Q, is 
net radiant energy exchange in the same units as the 
evaporation; is evaporation estimated by a simple 
equation of the Dalton type; 7, is temperature of the 
water surface; 7, is temperature of the air; ¢, is saturation 
vapor pressure at the temperature of the water surface; 
¢, is vapor pressure of the air above the water; and A is a 
factor applied to account for the difference between open 
water and a vegetated surface evaporation (Penman 
suggests 0.6 to 0.8 depending on the season; for the sake of 
simplicity, we use 0.7). 

The evaluation of the Penman formula, equation (1), 
involves estimating or measuring the energy received by a 
surface Because it is presumed, in the Penman approach, 
that the distribution of energy takes place according to the 
so-called Bowen ratio, no decision has to be made as to 
transfer of sensible heat to the water surface. 

Rijkoort [2] and van der Bijl [3] present nomographs for 
the solution of the Penman equation. This paper de- 
scribes a simpler graphical solution than van der Bijl’s 


homographs. 


2. GRAPHICAL SOLUTION 


The equation first derived by Penman [1] is used here 
with Units are that 
the formula will vield estimates of potential evapotrans- 
For the purpose of 


some modification. converted so 
piration in units of inches of water. 
these computations, the latent heat of vaporization is 
computed for approximately the average temperature 
observed during the summer in South Carolina. , in 
equation (1) is given by 


— §7T4(0.56—.092y e,)(0.10+0.90n/N) (2) 


where 2, is mean monthly extraterrestrial radiation in 
mm. day~! (evaporation equivalent); 7 is reflection coeffi- 
cient, used as 0.05 in Penman formula; n/N is ratio of 
duration of actual sunshine to maximum possible sun- 
shine; S is Boltzman’s constant =2.01 107° mm. 
(evaporation equivalent); 7, is air temperature in’ F.°; 
and ¢, is actual vapor pressure of air in mm. Hg. 

Fk, in equation (1) is estimated from a Dalton-type 
equation : 


=0.35(e,— e,)(0.5+ .0098 (3) 


where ¢, is saturation vapor pressure of the air in mm. 
Hg.; ¢, is actual vapor pressure of air in mm. Hg.; v% Is 
wind speed at 2 m. in miles per day; and £, is in units 
of mm. day7!. The constant 0.5 is used in lieu of the 
constant 1 that Penman used in 1948. 

By substitution of (2) and (3) into (1) and conversion 
of units so that the formula vields inches instead of 
millimeters of water, the Penman equation becomes [4 


(.7)(.0394) 


D937 DR,A—r)(.18+.55 n/N) 


(.7)(.0394) 
(.27)(.35)(e,— é,)(.5+ .0098 


D(.7)(.0394) 


ST,4(.56—.0924 €,)(.10+.90 n/N). 4) 
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Ficvre 1.—First of three graphs for computation of potential evapotranspiration by Penman formula. This graph yields F; as given 
by equation (6). To use the graph, (1) enter temperature on the upper left scale, (2) move horizontally to the radiation value or 
calendar date, (3) descend vertically to sunshine percentage, and (4) move horizontally and read value of FE; on the lower right 
stale. Dates shown on this graph are for Columbia, S.C, 
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Pictre 2. Second of three graphs for computation of potential evapotranspiration by Penman formula. This graph vields FE, as given by 
equation (7). To use this graph, enter temperature on the upper left scale, (2) move horizontally to the dew point (average), (3) 
descend vertically to the wind speed (daily average), and (4) move horizontally and read value of FE, on the lower right sesle 
Wind speeds for the Columbia, 8.C., anemometer (36-ft. height) are shown to the left of those for an anemometer at 2-m. height. 
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Figure 3.—Third of three graphs for computation of potential evapotranspiration by Penman formula. This graph yields FE as given by 
equation (8). To use the graph, (1) enter temperature on the upper left scale, (2) move horizontally to dew point, (3) descend verti- 
cally to sunshine percentage, and (4) move horizontally and read value of EF; on the lower right scale. Note that the values of Fy are 
negative. The computed value of potential evapotranspiration, E,, is given as the algebraic sum of F;, E2, and EF; obtained from 
figures 1, 2, and 3, respectively. 
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A graphical solution of the Penman formula can be 
devised from the three terms on the right side of (4). 
Thus 

E,=E,+E.+E, 
where 
(.7)(.0394) 


DR, —r)(.18+.55 n/N) 


(.7)(.0394) 


(.27)(.35)(e,—e,)(.5+ .0098 (7) 
) (-7)(.0394) 


+ .27 


E, €,)(.10+.90 n/N). 


(8) 


Figure 1 is the graphical solution of equation (6). 
Since 1) is a function of the air temperature, the upper 
part of figure 1 is drawn as temperature versus /?,, the 
extraterrestrial radiation, while the lower portion of the 
figure accounts for the influence of the variable n/N. 
The extraterrestrial radiation for any particular location 
is fixed by the date of the year. The reader may find 
complete tables of 2, values for all latitudes in the 
Smithsonian Meteorological Tables [{5|. It is possible, 
therefore, to adapt figure 1 to any fixed location by 
labeling the indicated isolines of 2, with the corresponding 
dates for that location. To use figure 1 the reader should 
enter the upper left hand margin with the air temperature, 
move to the right until the proper 7, is intersected, then 
descend to the n/N value. The solution F, is then found 
on the lower right hand margin. 

Figure 2 is the graphical solution of equation (7). The 
upper portion of this figure solves the actual and satura- 
tion vapor pressure functions. The relation at the bot- 
tom of the figure accounts for the effect of wind speed. 
To use figure 2, enter the upper left hand margin with 
the actual air temperature, move horizontally to the 
dewpoint, then descend to the wind speed curve. The 
solution /, is found on the extreme lower right hand 
margin. 

It should be noted that the wind speed wu, is that 
recorded at the 2-m. level. If this value is not known, a 
measurement at a height 4 may be converted by use of 


uz log h 
log 6.6 


log 6.6 
log h 


2 Up 
where the height of measurement is in feet. Thus this 
figure can be adapted to any particular site by relabeling 
the wind speeds with the corresponding values for the 
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particular anemometer height; e.g., for 36 ft. as applicable 
to the Columbia, 8.C., anemometer, given by the second 
set of wind speed labels. 

Figure 3 is the graphical solution of equation (8). The 
upper portion of this figure solves the temperature-vapor 
The lower part of the graph accounts 
lo use figure 3, 


pressure relation. 
for the influence of the parameter n/N. 
enter the upper left hand margin with the air temperature, 
move horizontally to the dew point, then descend vertically 
to the percentage of sunshine. The answer F is found on 
the lower right hand margin. 

The solution /, to equation (4) is then accomplished by 
adding results /,, ,, and /; found from figures 1, 2, and 3. 
For example, at Columbia, S.C., on May 1, with average 
temperature 80° F., dew point 70° F., wind speed 8 m.p.h., 
and sunshine 10 percent, figure 1 gives /,=0.072 in.; 
figure 2 gives E,=0.025 in.; and figure 3 gives /,= --0.011. 
Thus, from equation (5) 


in. 


3. CONCLUSION 


This paper presents a graphical solution of the Penman 
formula as adapted for Columbia, S.C. The graphical 
solution is given in the form of three graphs (figs. 1, 2, 3) 
which may be adapted for any location by changing 
figures 1 and 2 as described in the text. 
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upper-air winds has been checked by using available listings of vertical shear and extreme winds. 
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ABSTRACT 


The existence of undetected errors in recorded wind observations may have a biasing influence on a statistical 
In the progress of some studies it has been found necessary to reexamine the data being used. 


A series of 
The developed 


procedure permits correction for major errors and tolerates the minor (random) errors. 
The test of data by maximum wind profiles uses the highest and second highest scalar wind speed for each station 


and checks the data by profile sean. 


derived, exceedance of which marks the data as suspicious. 
In this program 3.5 percent of the observations ' proved suspicious and 85 percent 


suspicious value or it is corrected. 


thereof, that is, 2.9 percent of the observations, required correction. 


The test of data by vertical wind shear uses a critical value, theoretically 


A detailed check of the wind observation verifies this 


Thus the critical value is highly efficient. 


The errors were traced and split into clerical errors (1.1 percent), instrumental errors (1.3 percent), and computa- 


1. INTRODUCTION 


For use in missile design and performance studies by 
Army Ballistic Missile Ageney, basic upper-air wind 
observations were obtained for locations in the Pacific 
Ocean, North America, and Europe. The stations are 
listed in table 1. Preliminary analysis of the data revealed 
that there were occurrences of apparent errors in the 
observations as presented on punched cards. It was 
decided that these data should be checked. 

Although it was considered desirable to check and verify 
all the upper-air data, this was not possible because of 
the cost and the time required to review the mass of 
observations. Instead it was decided to establish a 
checking program which would permit a review eliminat- 
ing major errors, yet tolerating minor (random) errors. 
It was considered to be sufficient to restrict the checking 
process to higher wind magnitudes and wind shears, where 
the possibility exists that the reported extreme wind 
velocity arises from the addition of wind data and error 
with the same sign. 

This method permitted the correction of the major items 

“Formerly associated with the National Weather Records Center, U.S. Weather 
Bureau, Asheville, N.C. 


**Formerly associated with the U.S. Army Ballistic Missile Agency, Huntsville, Ala. 
Observation” used here refers to the entire ascent though sometimes only the value 


at a single level had to be corrected. 


tional errors (0.5 percent), which are quite within reasonable limits. 


(maximum wind speed and wind shear) required for 
missile design studies and verification of the wind data 
at the same time. 

The mass of punched cards was converted to magnetic 
tapes for use on high-speed electronic computers (IBM 704 
and 709). Use of these computers permitted the rapid 
searching of the data and machine listings of all observa- 
tions (plus the associated profile) which produced the 
higher wind speeds and wind shear for each altitude level. 
Also, it was possible to provide preliminary frequency 
distributions of the wind shear and speed data for use in 
further evaluation of suspicious data. The corrected 
observations were subsequently incorporated into the 
original data records and utilized on various statistical 
programs for use in missile design and employment 


studies. 


2. THEORETICAL BACKGROUND 


The problem existed that frequency distributions for 
wind shears and extreme values had been programmed, 
and tabulations similar to table 2 (described later) had 
been made before the necessity for critical review and 
correction became apparent. Thus, the problem was not 
to establish a suitable statistical theory of fitting extreme 
value data, but rather to develop an economical tech- 
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TABLE 1.—Station listing of verified wind data ! 


Job WBAN Period of Hour? Type of 
station Station name No. record (GMT) | equip- 
No ment 


i 


Cocoa (Patrick AFB) 4, 6, 


Cape Canaveral. 
Tokyo (Haneda AFB) 


Nagoya 
2411 
33123 


35140 
35010 


Fairbanks. . 
Tripoli (Wheelus AB) 


Berlin 

Wiesbaden 
34049 
13017 


Bitburg 
Port Lyautey 


Thule 17605 


56-12/57 
22536, 
12712 
14918 


Lihue 
Grand Bahama Isl 
International Falls 


Guam (Anderson AFB) 


San Juan. 11631 


5 | Balboa, C.Z. (Albrook AFB). 10701 | 


§ | Johnston Isl. AFB 21601 


Luzon (Clark AB 41207 
23236 
23273 
16204 
16201 


Santa Maria 
Keflavik (Meeks Fld). 


(Washington, 


Silver Hill 
D.C.) 


93722 


Adak 25704 


Fl Paso 23044 


Barrow 27502 
Alert 18601 


! The sequence of stations is not identical with that in tables 2 or 4. 
2 Prior to May 1957 observations were taken at 03 and 15 instead of 00 and 12 
not listed here. 
' Type of equipment: 2 
4 


This is 
Re abal, single theodolite 
C R-584 and SC R-545 
6 sc R-658 
s=GMD-1 and GMD-1A 


nique with minimum amount of time and cost using 
these available tabulations for selection and checking 
procedures. 

From the statistical point of view it must be pointed 
out that the problem of establishing profiles of maximum 
wind speed (magnitude of wind vector) or distributions 
of maximum wind shears should be approached otherwise 
than shown in this report by employing an extreme value 
statistic. However, frequency distributions of maximum 
wind speed and wind shear in the upper air are not com- 
pletely known and application of techniques similar to 
Thom’s [2, 3] method or Gumbel’s [1] theory would have 
necessary the careful evaluation of several such 
statistical This time-consuming basic study 
could not have been completed within the limitation of 
the available time and funds. 

The first part of the verification program, checking 
maximum wind profiles, was relatively simple. The listed 
and plotted profiles were scanned for suspicious values. 
Verification of those values considered suspicious was per- 


made 
systems. 
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formed by reference to the original data records. Details 
of the procedures employed are presented in section 3. 

The selection of suspicious values for wind shear data 
appeared very complicated in the beginning. Distribu- 
tions of wind component shear frequencies were available 
in the form of table 2 without the column marked ‘“Essen- 
wanger’s sum”. In the form shown in table 2, frequency 
distributions were given without regard to algebraic sign 
of the wind shear. The problem was to find a value 
which separated the acceptable values from suspicious 
values without reviewing too many observations or ac- 
cepting a large amount of unreliable data. This value 
will be called the critical value e,. 

The column marked 99.865” was practically identical 
with the maximum shear value. There is no reason to 
expect all maximum shear values to be wrong. On the 
other hand, the value listed in the column marked ‘97.72” 
was assumed to be acceptable since regularly it was ex- 
ceeded by 2 out of every 100 values. Thus we may 
contemplate the following ideas. 

We start with the assumption that zonal and meridional 
wind components are normal (Gaussian) distributions, or 
approximately normal. Departures from normality will 
be introduced later in the discussion. The wind shear 
data must then also follow a normal distribution. The 
shear distribution, disregarding the sign, then is a folded 
distribution. This folding occurs at the zero wind shear 
value. The mean value of the frequency distribution of 
zonal or meridional shear values, not disregarding the sign, 
usually will not coincide with zero. The question is now: 
Which portion of o (standard deviation) corresponds to 
the listed percentage frequencies, 50, 84.1, 97.72, ete.? 

If the folding occurs outside the +3 o value ? in refer- 
ence to the mean, then practically all shear values have 
the same sign. Then the shear distribution follows a 
normal distribution and the 50 percent value (median) 
virtually coincides with the mean. This statement holds 
for folding above +1.5 o (fig. la) if we assume that the 
frequency of data above 3 a is negligible. If the folding is 
within +1.5 o, then the mean value of the unfolded 
normal distribution must be smaller than the 50 percent 
value (fig. 1b). Thus the listed 50 percent value permits 
evaluation of the magnitude of this mean value. We recog- 
nize (see also table 2) that for practical purposes this value 
is so close to the zero wind shear that we can continue our 
discussion about the folded distribution as if it were 
folded at the mean value zero. Then the 84.1 percent 
value corresponds to 1.41 ¢, the 97.72 percent to 2.28 ¢, 
and the 99.865 percent to 3.20 co. 

We build the ratio 


3.200 
2.2860 


99.865 percentage value _ 


=1.40 
97.72 percentage value— 


Thus, theoretically we should expect the factor 1.4. A 


? Some statisticians may want to use the symbol s, the estimate of the population 
standard deviation «. This has no influence upon the development in this paragraph 
insofar as later the 84.1 or 97.72 value may express the empirical replacement for @ or s 
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10/52-11/56 00,1 
12868 | 11/56-12/57 00,1 
2 43311 | 1/51- 4/52 00, 1 2, 
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es | | 55-12/55 09, 2 : 
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T \ BLE 2.—Meridional wind shear (sec.-!) distributions at Keflavik, Iceland, for January. Station altitude, 54 m. MSL; latitude, 63°57’ N.; 
longitude, 22°37’ W.; period of observations, January 1951—December 1957 


| Number Cumulative Percentage Frequency | Essen- Suspi- 
Altitude km. (MSL) of OBS wanger’s _Maxshear | Pet. freq. cious 
15.9 50.0 99. 865 | sum . values 


. 0026 . 0090 ° ° . 036 . 0808 . 0621 
. 0008 . 0037 0291 0342 ‘ 
. 0030 . OO . 0314 
. 0001 . 0028 . 0176 0343 
. 0028 . 0306 
. 0037 . 0057 0319 
. 0004 . 0027 . 0163 
. 0003 . 0021 003 006 0267 
. 0004 0027 q 0724 | 
. 0023 006 ‘ O656 
. 0002 . 0025 . 004 : 57 0333 
. 0028 0480 
9.0- 10.0 . K . 0023 . 008 . 007 0345 | 
11.0-12.0__- 26: 0003 . 0383 
12.0-13.0__- 0277 
13.0-14.0 0256 
14.0-15.0_. 0242 
15.0-16.0 
16.0-17.0 
17.0-18.0 
18.0-19.0 


0046 
0027 
0073 
0008 
0074 
0008 


41.0 


review of the frequency tabulations showed that for the Psat +0.005 sec. 


meridional shear the average empirical value amounts to ‘ 
1.8 and for the zonal shear it is between 1.5 and 1.7. The where Pa, is the 84.1 percentage value and Pg;,72 the 97.72 
factor 1.6 seems, therefore, a sound compromise between Percentage value. 
theory and practice. This takes care of departures from 
the normal distribution law and a mean value different 
from zero. 

We have now established a theoretically acceptable 
critical value (e,) of 1.6 times 2.28 ¢, which equals 3.65 o. 
If the o is known, we can easily compute this critical 


value 

We also could use the 97.72 percentage value for the 
2.28 ¢ value. For the individual case, however, too many 
random variations may influence the result. Therefore, 
we may try to incorporate another procedure to decrease 
this effect. When we add the 84.1 percent and the 97.72 
percent values, we obtain 1.41¢ plus 2.280, which equals 
3.690. This is very close to 3.650. Thus we may 
derive the critical value by employing the 84.1 and 97.72 
percentage values. 

Further consideration may be given to an observation 
tolerance error. An error of 5 m. sec.~' per 1000 m. for 
those extreme values seems to be within the limitation of 
measurements.* Thus we tolerate this error for the 
critical value and derive finally oo tes 


FREQUENCY ([ Y } 


FREQUENCY Y 


Figure 1.—Folding of normal distribution and relation to 50 


Che first 3000 m. in table 2 are listed in 500-m. layers. It was decided for simplicity : ; 
percent line (median). 


to dept 0.005 sec.~! for those 500-m. layers, too. 


1 
\- 
n 
e 
e 
] 
D 
e 
21.0-22.0._. SY 0021 0042 0085 . 0176 0289 0311 0289 12 
ad 22.0-23.0... 75 . 0035 . 0064 0163 0269 0277 0269 33 | 
23.0-24.0 65 0017 . 0042 0074 0214 0290 0338 0290 
V 24.0-25.0 46 . 0019 0033 . 0058 . 0222 . 0255 0330 0255 17 ; 
25.0-26.0__- 36 0042 0092 0211 0211 0353 0211 78 
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The unit of «, is the same as in table 2, namely inverse 
seconds. From table 2 we may give a sample for compu- 
tation of the critical value for the layer, surface to 0.5 
km. It would be 


.0202 + .0369 .0050 = .0621 sec.~! 

This value is listed in table 2 in the column titled ‘‘Essen- 
wanger’s sum’. All shear observations exceeding this 
value as computed for each level were labeled suspicious. 
They are marked by an asterisk in table 2. Further 
description of the shear checking process will follow in 
section 4. Of the suspicious values, 85 percent had to be 
corrected, which is considered as high efficiency for this 
checking procedure. 

Results of randomly picked shear values for checking 
have demonstrated that efficiency drops sharply for values 
below the critical value. Thus, by the outlined process 
we have achieved the goal to eliminate major errors, and 
tolerate (random) minor errors which have little bearing 
upon the determination of missile design criteria. 


3. TEST OF DATA BY MAXIMUM WIND PROFILES 


The highest and second highest (scalar) wind speeds for 
each of the stations at each of the 45 levels (or to the 
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Maximum wind profile before verification, station no. 1. 
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highest level attained if it was less than 41 km.) were 
machine selected, and the observations * containing these 
high wind speed values were listed. 

The verification of these speeds may be managed in 
two ways: verify or correct each and every value, or locate 
and correct the greater majority of the erroneous values, 
particularly values that are very large or appear to be 
inconsistent with a smooth profile which would be ex- 
pected if all observations were correct. Practical economic 
considerations demanded that the latter be the guiding 
principle in the verification. One-third, or 8, of the sta- 
tions received total verification; i.e., every high wind 
speed value was checked to provide a basis of comparison 
to determine the adequacy of the smooth profile verifica- 
tion procedure. 

The verification procedure was divided into three steps 
or categories, namely: observation scan, terminating val- 
ues, and profile scan. The sequence of the checking 
process might have been arbitrarily established , however, 
it appeared that the sequence listed above would provide 
maximum assurance that the final product, the maximum 
wind profile, was correct. 


4 The word “observation” in this report is used in the sense of characterizing the entire 
ascent as one taken observation and is in this way different from an observed value. 
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Figure 3.—Maximum wind profile after verification, station no. 1. 
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Che observation scan was literally a visual scan of the 
observation. Each machine-selected observation was 
searched for apparent inconsistencies, such as a speed of 
10 m.p.s. followed at the next level by a speed of 116 
m.p.s., or similar rapid and large fluctuations in wind 
sped. Those observations containing such inconsistencies 
were checked and corrected as appropriate. This sub- 
jective scanning sufficed for errors detectable by discon- 
tinuity. Were all errors of this type, no further checking 
had been necessary. But though some profiles showed 
considerable smoothing from this process, most profiles 
still contained irregular contours which appeared sus- 
picious. The checking process was therefore continued. 
Each value occurring as the terminating speed for a 
given observation was considered suspicious. All observed 
wind values are the result of 2- or 4-minute averages 
except the last speed obtained which is allowed to be a 
l-minute calculation if a 2- or 4-minute average is not 
available. Terminating values were therefore checked for 
representativeness. This eliminated the error from termi- 
nating fluctuations, but further smoothing was needed. 
The profile scan necessitated the construction of the 
vertical profiles of these highest and second highest wind 
speeds. Figure 2 pictures the original profiles for station 
no. | to serve asan example. The first move was to apply 
to the profiles all corrections resulting from the first two 


MONTHLY WEATHER REVIEW 


201 


steps of the verification procedure. The profiles, after 
this preliminary correction, were then examined and ques- 
tionable values were “picked off” for verification. It is 
obvious that a great deal of subjectivity was also encount- 
ered here, but a system was utilized. The values were 
chosen in sets, each set being verified and the resulting 
corrections applied before the next set was chosen. Each 
set consisted of those wind speed values which, if changed, 
would most smooth the contours of the profile. This set 
checking process was continued until no errors were found 
or so few were found that the profile was virtually un- 
affected. 

Figure 3 shows the maximum wind profile after veri- 
fication. Since this investigation dealt with only the 
highest and next to highest speeds, and values to replace 
these when they were deleted or changed to a value below 
that of the second highest were not included in this pro- 
file, thus dashed areas of no data appear. This meant the 
highest or second highest wind speed for a dashed 
level would have to be obtained by going back to the 
original records and selecting the now highest (or second 
highest) value. As this would have had to be done by 
hand at National Weather Records Center, it was decided 
to leave the profile as in figure 3. The missing values 
may easily be replaced by machine selection at Army 
Ballistic Missile Agency. The corrected profile, in 
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general, takes on a smoother appearance with one distinct 
layer of maximum wind speed. This is true with all of the 
stations. 

Table 3 denotes the amount of checking and changing 
done in the maximum wind profile verification program. 
The data are subdivided into two groups; the first being 
for the 16 stations for which the verification was accom- 
plished by the procedure described in the preceding para- 
graphs, the second for the 8 stations which received total 
verification. The combination of these yields the overall 
results. Table 3 is self explanatory. 

As can be seen from table 3, there is some question 
concerning the efficiency of this three-step procedure, 
particularly when the percentage of values changed be- 
comes very large. It would seem that there still are a 
significant number of erroneous observations unchecked. 
It must be remembered, however, that the criterion for 
terminating verification was that the contours of the pro- 
file remain virtually unchanged, since the original goal of 
the verification was to obtain a representative maximum 
wind profile. This goal was attained in all cases where 
the profile was not completedly destroyed. Three sta- 
tions lacked sufficient data and for one the entire profile 
would have needed reconstruction. 

This method was designed for adequacy, with efficiency 
being second in importance. Methods aspiring toward 
high efficiency may be similar to the above procedures but 
would necessarily include more objectivity and less sub- 
jectivity. The method outlined for vertical wind shear 
verification with objective selection may be adaptable to 
wind speeds using similar frequency distributions. 


4. TEST OF DATA BY WIND SHEAR 


Selections were made from the shear distribution tables 
similar to table 2. These tables list the vertical shear by 
l-km. layers, beginning at 3 km. Below 3 km., 500-m. 
layers were taken. The first attempt at selection was by 
picking out values that seemed to be erratic. This was 
abandoned almost immediately since it was too laborious 
and uncertain (nearly 14,000 observations contained all 
of the maximum shear values). 

Then it was decided to use the method outlined in 
section 2. This method is as follows. Obtain the sum of 
the shear values at the cumulative percentage frequencies 
of 84.1 percent and 97.72 percent plus 0.0050 sec.~', 
then select shears for checking on the basis that any shear 
value greater than this sum is suspicious. In addition, a 
selection by observations was made: any value that had 
been obtained at a level reached ten or less times during 
the period of record was considered as suspicious, regard- 
less of magnitude. 

The selection of the suspicious values is demonstrated in 
table 2. This tabular form portrays layer versus cumu- 
lative percentage frequency with supplementary columns 
for observations counts and the observed maximum shear. 
The vertical wind shear (henceforth shear) values were 
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computed for both the zonal and meridional components. 
Listings of the profile for the maximum wind shear ac- 
companied the listings of the form of table 2. 

In the first set of tabulations maximum shear values 
only were subjected to verification. Later all shear val- 
ues exceeding the “critical” threshold (Essenwanger’s 
sum) were listed as a review program. It was understood 
that, regardless of what portion or portions of an obser- 
vation first raised suspicion, the entire observation was 
to be subjected to a checking process. 

Secondary review was limited to a review check program 
which was designed to encompass not only highly suspi- 
cious observations not contained in the first phase, but 
also observations in which errors may have occurred or 
been overlooked during previous verification. The latter, 
fortunately, occurred very few times. Even with this, 
some few errors undoubtedly still escaped detection. 

The combination of these methods has proven accept- 
ably efficient in that 85 percent of the observations 
selected as suspicious were found to be erroneous. 

The errors have been tabulated in three categories: 
clerical, instrumental, and computational. The clerical 
errors were based on the premise that no technical ob- 
server training was required to perform the work classified 
as clerical. They were subdivided into three types: 
punching, extraction (or transcription), and plotting of 
Form WBAN 20A. The punching was, of course, the 
production of the cards constituting the original card 
decks. The extraction was the “picking off” of values 
from the WBAN 20A, a plot of the observed data in the 
form of wind speed and direction versus height. Plotting 
refers to wrong plots of WBAN 20A. 

Instrumental error may be defined as large and rapid 
fluctuations in angles (azimuth and elevation) incompat- 
ible with the calculated height changes. 

The computational error was subdivided into two types: 
Calculations on WBAN 20 (the observer’s work sheet) 
and fictitious ascension rates (of the balloon) as calculated 
from erroneous pressure-temperature-time measurements. 
Fictitious ascension rates may be thought of as instru- 
mental, but only the tracking equipment was considered 
an instrument in this study. The tracking equipment 
consisted mainly of the theodolite (visual tracking), the 
SCR-658 (manual radiosonde tracking), and the GMD-1 
and GMD-1A (automatic radiosonde tracking). 

Table 4 shows the error statistics resulting from the 
verification of the shear selected observations. The head- 
ings are self explanatory. Some observations contained 
more than one type of error so that the occurrence of 
errors exceeds the number of observations changed. The 
end results in this tabulation prove interesting in the 
predominance of the instrumental and clerical errors. 
The relatively small computational error is gratifying. 

We note in table 4 that from all observations only 3.4 
percent were found to contain suspicious values; from the 


suspicious observations 85 percent, i.e., 2.9 percent of all 
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« Sequence of stations is not identical with station listing in table 1, 


observations, had to be changed. This means that one 
or several values had to be corrected in 2.9 percent of the 
observations. As one observation contains numerous 
levels, the actual percentage of errors in relation to the 
wind data of all levels is far less. The clerical error con- 
tributed 1.1 percent, the instrumental error 1.5 percent, 
A perfect card 


As one should 


and the computational error 0.5 percent. 
deck would be the ideal goal, of course. 
expect, this goal cannot be reached without a thorough 
check of the produced data. Thus, the magnitude of 
percentage of the larger error discovered by the method 
outlined above stavs well within reasonable tolerance 
limits. 

A few brief remarks may be in order concerning the 
First, the major 


division of errors is not as clear-cut as it appears. Prior 


difficulties involved in’ verification. 
to 1956, the extracted values were not always entered 
in the allotted space on WBAN 20A. Errors occurring 
subsequent to plotting were attributed to extraction in 
such cases. Then, beginning on January 1, 1956, WBAN 
20\ was no longer a filed form so that plotting errors 
were undetectable. WBAN 20B had come into being, 
on which was allotted space for the ‘“‘extracted”’ values. 
Again, arbitrarily, errors occurring between calculation 
WBAN 20) and punching were attributed to extraction, 
80 that an overlap occurred. Fortunately, plotting errors 
contributed very little to the mass of clerical errors. 
Perhaps the most difficult problems arose with erratic 
instrument data and fictitious ascension rates. Balloon 


determined through pressure-temperature 


heizhts are 


considerations on Forms WBAN 31 A and B. Errone- 
ously low pressure or high temperature readings result in 
determinations, which when 


present unduly 


high height 


with 


erroneously 
ascension time, 


This normally 


considered ay 


large ascension rates. occurred in the 
topmost lavers of the higher ascents. This same reason- 
ing can also account for highly erratic or incompatible 
height changes at any level with the erroneous pressure 
or temperatures high or Highly 
angles are caused mainly by equipment malfunction and 
limitation. Any turbulence 
would be confined to the lowest levels except under rare 


being low. erratic 


appreciable influence — of 
conditions such as balloon entry into a thunderstorm or 
possibly clear air The first event 
eliminated by consideration of the weather reports. 


was 
The 


latter event, generally not too frequent in occurrence, is 


turbulence. 


not known at the present in sufficient detail to make an 
unequivocal decision. It was felt, however, that when 
low ascension rates are involved, it is more likely that 
the data are erroneous due to instrumental errors. This 
conclusion may be due for revision after knowledge of 
clear air turbulence has improved. 

In general, erratic data were determined subjectively 
since there were no adequate objective methods available. 
Ascension rates were, in general, declared fictitious in a 
subjective manner since time limitations and unavailable 
observation data precluded thorough checking. There 
are two compensating factors, however, in that observa- 
tion data were always given the benefit of the doubt, and 
all decisions were made by qualified meteorologists. 


61 Jur 1961 
its, 
1@s 1 2 3 5 6 7 
rs 
er- 
1 5,114 159 167 42 57 
2 4, 992 | 556 223 104 249 101 93 ‘ 
3 5,114 | 537 169 133 134 27 | 100 | ' 
4 5,114 | 26 208 159 ri.) 58 
5 3, 652 5M) 106 87 107 33 | 
pi- 
ut 
or 
er, 
is, 
ns 5,114 127 116 121 lod 
as 21 5,114 48h 217 108 227 47 139 
% 5,114 521 12 146 164 30 113 
pg: 4 3, 652 502 103 93 35 49 
{ 
al 
ed 
PS: 
of 
he i 
rd 
les 
he 
ng 
vid 
it- 
38: 
ft 
ed 
ts. 
u- 
ed 
nt 
he 
d- 
ad 
of 
he 
he 


204 
5. CONCLUSION 


This study has discussed the possibility of checking 
wind data by maximum wind profiles and wind shear dis- 
tributions. While the maximum wind profiles were evalu- 
ated for suspicious values by profile sean, the checking 
process by wind shears was based upon computation of a 
critical value e.. Exceedance of this critical value made 
the shear value suspicious and subject to verification. 

The derivation of the critical value «, was developed 
and the application to 24 stations showed an efficiency of 
85 percent, which may be considered very high. Although 
3.5 percent of the observations proved to be suspicious, 
and 2.9 percent had to be corrected, the actual corrections 
are less, as one observation in the average contains be- 
tween 20 and 30 level values, not all of which had to be 
corrected, 

The errors were traced and divided into 1.1 percent 
clerical errors, 1.3 percent instrumental errors, and 0.5 


percent computational errors. These are within reasonable 


limitation, 

It may that 
wind profiles or maximum wind shears may be better ap- 
proached by theoretical statistical processes in order to 


be stressed establishment of maximum 
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eliminate the effect of the relatively short period of avail- 
able data record. Time and cost limitations, however, 


prevented further investigation in this direction. 
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1. HIGHLIGHTS 


The mild weather of the second half of February 1961 
[1] persisted into March as temperatures averaged above 
seasonal normals in nearly the entire United States (fig. 1). 
A notable exception was Alaska, where near record cold 
prevailed. Frequent storminess in the central Plains, 
associated with a mean trough (fig. 2), brought heavy 
precipitation and flooding to much of that area. This 
March was in sharp contrast to March of 1960 [2], a very 
cold, snowy period in the eastern two-thirds of the Nation. 


2. AVERAGE CIRCULATION 


The monthly mean circulation patterns at 700 mb. 
(fig. 2) and sea level (fig. 3) for March 1961 were featured 
by strong blocking in the Atlantic. This was associated 
with a band of positive anomaly of both 700-mb. height 
and sea level pressure around 55° N., and a much deeper 
than normal trough near mid-Atlantic. In addition, the 
Icelandic Low was split and displaced southward, a feature 
commonly associated with blocking. 

Blocking was also present over North America, where 
700-mb. heights averaged as much as 160 feet above normal 
over James Bay (fig. 2). A deeper than normal trough 
extended from mid-United States southwestward. The 
rather short wave spacing that existed at middle and low 
latitudes across North America and the Atlantic was quite 
characteristic of blocking. 

Blocking also appeared in the Pacific, where it was 
related to the center of positive height anomaly near 
Kamchatka, the deep trough in the Gulf of Alaska at 700 
mb. (fig. 2), and the southward displacement of the split 
Aleutian Low at sea level (fig. 3). 

The circulation in polar regions was dominated by a 
strong High centered north of Alaska at sea level and 
aloft (figs. 2 and 3). This cell developed strongly near 
mid-month, concomitant with appearance of a deep 
evclonie center in northern Scandinavia. Both cells re- 
mained quasi-stationary, reaching their greatest strength 
during the 5-day period March 23-27 when the maximum 
mean sea level pressure in the High was 1047 mb. and the 
minimum in the Low 982 mb. 


! Descriptions of the weather of April, May, and June 1961 will appear in the July, 
Auvust, and September issues of the Monthly Weather Review, respectively. 
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Over Eurasia the mean circulation was well defined, con- 
sisting of a deep trough in eastern Europe flanked by 
strong ridges, and there was little displacement of these 


features from sea level (fig. 3) to 700 mb. (fig. 2). The 
European ridge-trough system was associated with the 
strongest height anomalies in the Northern Hemisphere 
(+390 ft. and —470 ft. in fig. 2). 

The monthly mean isotach chart and its departure from 
normal (fig. 4) further reflect the blocking character of the 
circulation at 700 mb. The jet axis at this level was well 
defined and displaced south of its normal position over 
nearly the entire western portion of the Northern Hemis- 
Note also the unusual double jet across the At- 


phere. 
These two jet axes merged in the northeastern 


lantie. 
Atlantic to produce an extensive area of above normal 
wind speeds, as much as 11 m.p.s. above normal over 
southern Scandinavia (fig. 4B). 

Wind speeds well above normal were also observed in the 
eastern Pacific (fig. 4B). As this strong current ap- 
proached the west coast of the United States, it underwent 
marked diffluence, aided in part by blocking in Canada. 


3. CIRCULATION TRANSITION WITHIN THE MONTH 


The most important change in circulation during the 
month was associated with retrogression of a pronounced 
blocking wave which developed from a strong ridge over 


“4 = 
Figure 1.—Departure of average temperature from normal (° F.) 


for March 1961 (from [5]). Mild weather dominated the Nation. 
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FiGurReE 2. 


Mean 700-mb. contours (solid) and departures from normal (dotted), both in tens of feet, for March 1961. 


Blocking at middle 


latitudes was prominent over North America and the Atlantic, 


Europe in February [1]. During Mareh the blocking 
spread across the Atlantic into North America, accom- 
panied by an abrupt transition from high to low zonal 
index conditions (fig. 5). The magnitude of this transi- 
tion in circulation is shown in figure 6, which gives the 
change in 700-mb. height between the first and last halves 
of the month. Note the pronounced increase in heights 
at higher latitudes and falls at lower latitudes over the 
Atlantic and North America. The greatest change oc- 
curred in the Atlantic where 700-mb. heights rose as much 
as 740 feet. 


A further manifestation of this reversal in circulation 
is obtained by comparing typical 5-day mean patterns of 
700-mnb. height and sea level pressure observed at time of 
high index (fig. 7) and low index (fig. 8). Of additional 
interest are the tracks of primary, migratory eyclones 
(dashed) and antievclones (dotted) appearing during 


these 5-day periods and superimposed on the sea level 


charts. The high-index patterns (fig. 7) display most of 
the characteristic features of this type of circulation '3]; 
i.e., fast zonal westerlies, small wave amplitude aloft, and 
pressure systems oriented east-west. 
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FicgurRE 3.— Mean level isobars (solid, with intermediate 
isobars dashed) and departures from normal (dotted), both in 
millibars, for March 1961. Unusually strong pressure gradi- 


ents prevailed across Europe and Alaska. 


sea 


During the period of lowest index the circulation was 
characterized by a disruption of the westerlies into closed 
cellular centers, shorter wavelengths aloft, and a general 
north-south orientation of pressure systems (fig. 8). In 
addition, the split storm track in the Atlantie during low 
index contrasts with the single primary track across the 
Atlantic during high index. Note also (fig. 8B) the high 
frequency of slow-moving storm centers in the central 
United States with low-index conditions. Strong block- 
ing in the Atlantic was evidenced by development of an 
upper-level “omega” pattern (fig. SA). Retrogression of 
blocking into North America was accompanied by devel- 
opment of a similar “omega”’ pattern in the United States, 


although not as pronounced as in the Atlantic. 

An additional comparison of these extreme circulation 
states can be made by reference to figure 9, which shows 
the 700-mb. zonal wind speed profiles for the two periods. 
The sharply peaked profile associated with high index is 
in marked contrast to the double peaked profile and dis- 
placed westerlies with low index. 


4. WEATHER IN THE UNITED STATES 
TEMPERATURE 


Temperatures in the United States during March 1961 
averaged warmer than normal over all but the extreme 
Northeast, the interior valleys of the Pacific Coast States, 
and Alaska (fig. 1). Greatest positive departures were 
observed in the Northern Plains and Upper Mississippi 
Valley—as much as +10° F. in the eastern portions of 
Montana and North Dakota. Persistence of the tem- 
perature pattern, usually considered in terms of the total 
zero-plus-one-class change [4], was extremely high from 
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Figure 4.—(A) Mean 700-mb. isotachs and (B) departures from 
monthly normal wind speeds, both in meters per second, for March 
1961. Solid arrows in (A) indicate principal axes of maximum 
wind speeds, and dashed arrows their normal March positions. 
Regions with wind speeds greater than 16 m.p.s. and anomalies 
greater than 8 m.p.s. are hatched. The westerlies were south of 

their normal position over much of the hemisphere. 


February to March, as 93 percent of the country did not 
change by more than one class (out of five). This com- 
pares with an average of 68 percent for the period 1942-60. 

The mild temperature regime was quite well related to 
the mean circulation patterns. At 700 mb. the anomalous 
flow between the east coast ridge and the southwestern 
trough was predominantly southeasterly over all but the 
Far Southwest and Northeast, where northerly flow was 
associated with cooler conditions (fig. 2). Moreover, 
strong Pacific westerlies (fig. 4) combined with a deep 
trough in the Gulf of Alaska and southerly anomalous 
flow in western Canada (fig. 2) to result in relatively mild 
Pacific air masses throughout most of the Nation. The 
thickness in the layer from 1000 to 700 mb. also corre- 
sponded well with the temperature pattern (figs. 1, 10). 


96) 1961 
>). 
Pity 
\ aK 
- Ke Pps; 
if 
of 
il 
. 
| 
if 
d 


MONTHLY WEATHER REVIEW 


ro) 


METERS PER SECOND 
@ 


4 
23 16 23 30 6 


FEB MARCH 1/961 APR 


PiGuRE 5. 

over the western half of the Northern Hemisphere between lati- 
and 55° N, 
index values (plotted at middle of period and computed thrice 
weekly 
averaged from maps of [6]. 


March. 


Time variation of speed of 700-mb. westerlies averaged 


tudes 35 Solid line connects 5-day mean zonal 


while dashed line gives the corresponding “normal” 
Note pronounced index evele during 
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Figure 6.—-Change of mean 700-mb. height (tens of feet) from first 
to second half of March 1961. 
and falls at lower latitudes over the Atlantic and North America 


Height rises at higher latitudes 


reflect retrogression of blocking. 
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Figure 7.—Five-day mean charts for February 28-March 4, 1961 
for (A) 700-mb. contours (tens of feet) and (B) sea level isobars 
(millibars). Tracks of important migratory cyclones (dashed 

and anticyclones (dotted), with March 2 positions indicated by 

large dots, are shown on (B). These circulation patterns were 
associated with high-index conditions early in the month, & 


A good example of the effect of cloudiness and precipi- 
tation on the temperature pattern was observed during 
March in the Plains States. 
figures 1 and 11. 
sunshine were related to the large positive temperature 
departures in the Northern Plains, while to the south heavy 
precipitation, much of which was in the form of frequent 


This is seen by comparing 
Unusually dry conditions and abundant 


snows, Was associated with temperatures close to normal. 
The most unusual warmth occurred during the period 


March 4 to 6 in the eastern half of the Nation, when 
early-season maximum temperature records were either 
established or equaled at Peoria, Ill., Augusta, Ga., 
Wilmington, N.C., and Richmond, Va. This warimth 
occurred when the mean ridge over the east coast (fig. 2) 
reached its greatest strength. 
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Figure 8.—Five-day mean charts, (A) 700 mb. and (B) sea level, 
for the period March 18-22, 1961, with March 20 positions of 
evclones and anticyclones shown by large dots. These circulation 


patterns were characteristic of low-index conditions late in the 


month, 


Near the end of the month retrogression of blocking 
into Canada, along with below normal upper-level heights 
inmost of the United States, was accompanied by a gradual 
change from a warm to a cold regime. 

While the contiguous United States had mild weather, 
Alaska experienced an unusually cold March in all but the 
southeastern coastal area (fig. 1). A partial explanation for 
this extreme cold lies in the stronger than normal north- 
easterly flow from the Arctic Basin at sea level (fig. 3), 
along with below normal 700-mb. heights (fig. 2). In 
addition, the thiekness in the layer from 1000 to 700 mb. 
(fig. 10) was well below normal. Note also the unusual 
extent of the area of cold thickness at higher latitudes, 
With a center of extreme cold over the Davis Strait and a 


strong gradient to the south. 
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Figtre 9.—Five-day mean wind speed profiles for high- and low- 
index stages of the index cycle during March 1961. 
speed of the westerlies at middle latitudes was compensated for 


by an increase at lower latitudes, thereby tending to conserve the 


Decrease in 


total westerly momentum, 
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Figure 10.— Mean 1,000-700 mb. thickness departures from normal 
(tens of feet) for March 1961, with areas of subnormal values 
shaded. Abnormally cold air covered the higher latitudes of 

North America, but warm conditions were dominant elsewhere 


over the continent. 


PRECIPITATION 


Precipitation over much of the Nation was adequate 
during March, and more than twice the normal amount 
fell in portions of the Pacifie Northwest and the Central 
Plains (fig. 11). At Yakima, Wash., this was the second 
wettest March of record, precipitation being more than 
Near record amounts also fell 
This heavy pre- 


four times the normal. 
elsewhere in Washington and Oregon. 
cipitation was related to frequent storms moving eastward 
from the Pacific in faster than normal southwesterly 


flow (figs. 2, 3, 4). 
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Figure 11. 
(from [5]). 


Percentage of normal precipitation for March 1961 


Central Plains and Pacific Northwest. 


Heavy amounts also fell in the middle of the Nation 
where Des Moines, lowa and Kansas City, Mo. both re- 
ported record precipitation totals for any March. This 
resulted in considerable flooding in Iowa and portions of 
Minnesota and Wisconsin. Record high river stages were 
reached late in the month along the Cedar River in Iowa 
at Charles City and Waterloo. This heavy precipitation 
was associated with the deep mean trough in the South- 
west and related storm systems which moved from the 
southern Plateau to the Great Lakes region. 

Unusually dry conditions prevailed in the Northern 
Plains (fig. 11), where near record dryness was observed 
in North Dakota. 


storm 


This area lay between two primary 
tracks, one to the south 
Ridge conditions and northwesterly 
flow aloft (fig. 2) also were related to the dryness. 


and the other across 


southern Canada. 


Near record amounts of precipitation fell in the 
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Precipitation was also well below normal in souttierp 
portions of California, Arizona, and Texas (fig. 11). At 
Prescott, Ariz., the period from November 1960 to Mareh 
24, 1961 was the driest such period on record since 1865, 
Near the end of the month, however, a record 24-hour 
snowfall of 9 inches broke the long dry spell. The drvnegs 
in California and Arizona was related to stronger than 
normal northerly flow in the rear of the mean trough 
(fig. 2). The precipitation deficit in southern Texas is 
more difficult to relate to the mean circulation patterns 
at sea level (fig. 3) and 700 mb. (fig. 2). However, a 
possible explanation lies in the position of the 700-mb, 
mean jet maximum through central Texas (fig. 4A), with 
relatively heavy amounts of precipitation observed north 
of the jet and light amounts in subsidence and anti- 
cyclonic vorticity to its south. 
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